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Abstract 

We consider an application involving a financial quadratic portfolio of options, when the joint under- 
lying log-returns changes with multivariate elliptic distribution. This motivates the need for methods 
for the approximation of multiple integrals over hyperboloids. A transformation is used to reduce the 
hyperboloid integrals to a product of two radial integrals and two spherical surface integrals. Numerical 
approximation methods for the transformed integrals are constructed. The application of these methods 
is demonstrated using some financial applications examples. 

1 Introduction 

Value-at-Risk (VaR) is considered to be one of the standard measures of market risk. VaR measures the 
maximum loss that a portfolio can experience with a certain probability over a certain horizon, for example, 
one day. Mathematically, if the profit or loss is given by II(t) — IT(O), for which is the price of the 
portfolio at t, VaR for a confidence level 1-a is determined by the following equation: 

Pro6{IT(0, 5(0)) - U(t, S(t)) > VaR a } = a, 

where S(t) = (Si(t), . . . , S n (t)) is a vector of asset prices that govern risk factors. 

In this paper, we consider numerical methods for the estimation of integrals over hyperboloids and 
their application to VaR computations for a complex portfolio that contains options depending on the 
market fluctuations that create risk. We reduce the problem of VaR computations to multiple integrals over 
hyperboloids, and show how these integrals can be approximated using techniques described by Genz and 
Monahan 7 and Sheil and O 'Muircheartaigh|l 1| . 

One of the most important analytic methods for VaR computation, which is called A-normal VaR, was 
introduced in the RiskMetrics Technical Document (1996). The method is based on the assumptions that 
the distribution is Normal and the portfolio is linear. Sadefo-Kamdem j^j, generalized the A-Normal VaR 
by introducing the A-elliptic VaR for a linear portfolio, with the A-Student VaR given as an example. An 
advantage of the A-clliptic VaR (for example, A-Normal or A-Student) is that the formula is still fairly 

•Paper based on the M. Jules SADEFO KAMDEM University of Reims Phd Thesis. 



1 



simple to calculate. But in practice, if we deal with a A-hedged portfolio, the A-elliptic VaR does not 
provide a realistic model, and that is why alternatives were proposed in the paper of Brummelhuis, Cordoba, 
Seco and Quintanilla [2]. In that paper, mathematical stationary phase techniques were used to estimate 
VaR for a quadratic portfolio when the risk factor changes with a Normal distribution. In a sequel paper, 
Brummelhuis and Sadefo-Kamdem |3] provided an estimation method with more precision for a VaR with a 
quadratic portfolio and generalized Laplace distributions. 

We assume, as in 0, that the approximation for the price of the portfolio is given by 



and we also assume that the joint log-returns X is elliptically distributed. For further details about elliptic 
distributions, see Embrechts, McNeil and Straussman [5]. T, O and A are functions of some sensitivities of the 
portfolio (see Taleb [Tl|, 1997, for a discussion concerning sensitivities). We also suppose that t = 1, because 
the time horizon for VaR is generally taken to be one day If the log-returns of X are Normally distributed, 
Albanese and Seco [I] have shown how to reduce the analysis of a quadratic VaR to the computation of 
the integral of a Gaussian over a quadric in a space of possibly very high dimension. We will use the more 
general assumption of an elliptic distribution for the risk factors. 

This paper proposes numerical methods for the approximation of integrals over hyperboloids, with ap- 
plication to estimate the VaR. We combine some techniques described in [H] ,0 and JI] to approximate the 
integrals over hyperboloids for VaR, with the generalized assumption that the underlying joint log-returns 
changes with an elliptic distribution. To illustrate our method, we will take the familiar case of Normal 
distributions, and we consider test examples for two A-hedgcd portfolios from the French CAC 40 market. 
Brummelhuis, Cordoba, Quintanilla and Seco |2J have considered a quadratic portfolio with an analytic 
approximation for the Gaussian integrals over quadrics. Sadefo-Kamdem and Brummelhuis 3 have pro- 
vided a similar analysis with a Generalized Laplace distribution. Albanese and Seco pQ investigated the 
approximation of integral of a Gaussian over a hyperboloid region with Fourier transform methods. 

One of the most common methods for quadratic perturbations of the linear VaR uses the Cornish-Fisher 
expansion for the quantile function of non-Gaussian variables. There are also some quadratic approximations 
in Hull [B] and Dowd 0j. Many papers in literature have proposed numerical methods for the quadratic 
approximation (see, for example, |12j . where Sadefo-Kamdem proposed the use of some numerical methods 
of Genz, jjjj, and the use of hypergeometric functions for a portfolios of equities VaR with multivariate 
t-Student distribution) . 

The rest of our paper is organized as follows. In Section 2 and 3, following Albanese and Seco P, we 
show how portfolio volatility can be used to reduce the calculation of VaR to the approximation of integral 
over hyperboloid, assuming elliptic distributions that admit a density function. In section 4, we propose 
a numerical method for the approximation of integrals over hyperboloids using some methods of Genz UJ, 
Genz and Monahan [7] and To illustrate our method we use examples where the density function is 

Normal. In section 5, we consider two examples of financial portfolios, and we have showed that our method 
is applicable to estimate the VaR for the portfolio. In Section 6 , we provide some conclusions. 



2 Quadratic Portfolio of Options Application 

In this section, we will define a quadratic portfolio of options as Quintanilla did in ^U]. We first define 
X = (Xi, . . . , X„), with 

X, = log(Si(t)/Si(0)), 

and we define Ai = (A}, . . . , A"), with Aj = 5i(0).A* , and A = (A 1 , . . . , A"), the gradient vector of the 



n(t) - n(o) w9f + ax* 





n 



by 
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with T = (t i A = ( g f 2 "g. (0)) , the Hessian of the portfolio at time t = 0. If we use a 2 nd 

order Taylor series approximation for II , then 

u(t, s(t)) - n(o, s(o)) w te + AiX 4 + xriX*/2, 

where 6 = ^(0) . 

If we consider a A-hedged Portfolio, we have A = 0, and therefore Ai = 0. Our goal is to determinate 
the Value-at-Risk quantity V with confidence level 1 — a when t = 1, as a solution to the equation 

Gt, (-V) = P(6 + XTiX 4 /2 < -V) = a. 

If we assume that the joint underlying log-returns X have a multivariate elliptic distribution with zero 
mean, then Gti(— V) is given by 

dx 

' {Q+xr 1 x t /2<-v} \/ det(E) 



G Tl (-V)= f g{xY,- 1 x t )^^==a. 
Jit 



3 Transformation to a Hyperboloid Integration Region 

We first decompose E as £ = CC*, where C is the (lower-triangular) Cholesky decomposition factor of S, 
and then we use the transformation x = yC* to give 

Gr(-V) = [ g{yy t )dy. 

J {e+yCTCy* /2<-V} 

We next assume that the sensitivity-adjusted variance-covariance matrix, C*TC, has a diagonalization in 
the form C*TC = ODO*, with O orthogonal, and D diagonal. CTC* is not necessarily definite, but we can 
assume that a diagonalization has been constructed with 



where 

D e = 



D+ 
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the transformation z = yO to give 



< d+ < . . . < d+ . Then, we can use 



Gr(-V) = f g(zz*)dz, 

J {Q+zBz* /2<-V} 

and finally, we can use the transformation w — IDI 1 / 2 ^ to give 

Gr(-V)=[ g (HB|-V) (1) 

where w = (w + , w_ ) is the decomposition of E™ into the respective positive and negative subspaces of the 
eigenbasis for C'rC. After changing the direction of the inequality, we obtain the following expression for 
G(R), which will be the starting point for the discussion of our computational methods. 

G(R)=I 5 HB|-V)^=, (2) 

where R 2 = 2(V + &). The integration region is the hyperboloid defined by |to_ | 2 — \w + \ 2 > R 2 . Our goal is to 
determine R a solution to G(R) = a. Once we find R, we will have the approximate quadratic Value-at-Risk 
given by V = R 2 /2 - 6. 



3 



4 Integration over Hyperboloids 

4.1 The Normal Case 

For many applications, the distribution g is a Normal distribution. In these cases, 

dw 



G(R) = I 



e 



(|™-P-h + | 2 >;? 2 } y/(2n) n det{\B\) 
Separating the w variables, we find 

dw- dwj 



G(R) = / e- w + D + w +' 2 I e 



{|u)+| 2 >o} J{\w-\ 2 >B?+\w + \ 2 } \J (2n) n -det(D_) yj '(27r) n + det(D + ) 

The inner integral for u>_ can be efficiently computed using the algorithm described by Sheil and O'Muirchcartaigh 
[TT] . so we define H(R,r) by 

H(R,r)=[ e - W -n-^-,2 dw. 



'{|™_| 2 >i? 2 +r 2 } ^(27zyidet(D_) 
i 

and let w+ = D^_z Then can rewritten as 

/CO pOO pOC 
/ .../ e-^H&lzD+z*])- 
-OO «/ — OO J — OQ 



dz 



Integrals in this form can be approximated using methods described by Genz and Monahan [Jj. 
4.2 The General Case 

We need to determine approximations to integrals in the form 



G(R) = / ip(x,y)dydz, 
J{\*\ 2 ~\y\ 2 >R 2 } 

where x € R™ 1 and y € l" 2 . If we use the changes of variable: y = ^£,2, x = ri£i, with r 2 — \y\ and T\ = |x| 
then G(R) becomes 

/•OO /• />oo P 

G(R)= r^- 1 / r? 1 " 1 / v{ri^r 2 £,2)da^ 1 )dr 1 da{£ 2 )dr 2 



J y/W+T* J\d\=i 

We now have G{R) defined in terms of a product of two integrals over hyper-spherical surfaces, defined by 
(£2! = 1 and |£i| = 1, and two radial integrals. The hyper-sphere surface integrals can be approximated 
using methods described in the paper by Genz [S]. If the surface and radial integrals are combined, then 
generalizations of the methods described by Genz and Monahan [7] can be used. Efficient approximation of 
the radial integrals will depend on information about the rate of decrease of the integrand ip for large values 
of n and T%. 

5 Application Examples 

We will distinguish 3 case in our analysis : 
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• n_=0; if q is Normal, 



G(R) 



dw+ 



{\w+\*<R?} y/(2n) n +det(D + y 



G(R) can be efficiently computed using the Sheil and O'Muircheartaigh JT] algorithm, when R 2 = 
-2(V + 6)>0. 



n + =0; if g is Normal, 



G(R) 



-w-DZ wt/2 . 



dw- 



i{\w-\*>B?} y/(2Tr) n - det(D_) 

G(R) can be efficiently computed using the Sheil and O'Muircheartaigh 11 algorithm. 
• n_ and n + are both nonzero; g is Normal, 



/oo i>oc roc 
/ .../ e- x * t ' 2 H(R,\zD + z t \) 
-oo « — oo J —oo 



dz 



Integrals in this form can be approximated using methods that are a combination of the Sheil and 
O'Muircheartaigh Ilj and Genz and Monahan [Jj algorithms. 

5.1 An Example when n + =0 

We construct a A-hedged portfolio that contains n equities and n European call options on these equities 
from the French CAC 40 Market. The Price of the Portfolio is given by 

n 

U(t, S(t)) = J2[-Ci(t,Si{t)) + A* • Si], 

i=l 

where Si is an equity price i, with S(t) — (Si, . . . , S n ), and Ci(t, Si(t)) is the price of European call option 
i on equity i. A is known in the literature as a gradient portfolio sensitivity vector. Our portfolio has been 
chosen so that A = 0, with A* = §§ i (»S'i(0)), and A = (A 1 , . . . , A"). The exercise price of each European 
call option is given in the following table for an example where n — 9.: 

Table 1: Data for Nine CAC 40 European Call Options 





Exercise Price 


Interest Rate 


Maturity 


Underlying Price 


BNPPARIBAS 


44.26 


0.1 


3 month 


39.75 


BOUYGUES 


23.49 


0.1 


3 month 


27.30 


CAP GEMINI 


34.71 


0.1 


3 month 


24.00 


CREDIT AGRICOLE 


17.36 


0.1 


3 month 


14.80 


DEXIA 


11.5 


0.1 


3 month 


9.38 


LOREALL 


61.85 


0.1 


3 month 


62.90 


TF1 


26.38 


0.1 


3 month 


22.02 


THOMSON 


15.22 


0.1 


3 month 


17.13 


VIVENDI 


16.19 


0.1 


3 month 


17.00 
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Using the above data with the exponential moving weighted average (EMWA), we obtained the following X: 



V 



0.0017 
-0.0001 
0.0012 
0.0005 
0.0008 
0.0008 
0.0008 
0.0002 
0.0002 



-0.0001 
0.0009 
0.0005 

-0.0001 
0.0000 
0.0000 
0.0000 
0.0006 
0.0005 



0.0012 
0.0005 
0.0038 
0.0006 
0.0011 
0.0008 
0.0014 
0.0006 
0.0007 



0.0005 
-0.0001 
0.0006 
0.0006 
0.0002 
0.0004 
0.0004 
0.0001 
0.0001 



0.0008 
0.0000 
0.0011 
0.0002 
0.0015 
0.0007 
0.0008 
0.0000 
0.0002 



0.0008 
0.0000 
0.0008 
0.0004 
0.0007 
0.0011 
0.0006 
0.0004 
0.0005 



0.0008 
0.0000 
0.0014 
0.0004 
0.0008 
0.0006 
0.0013 
0.0000 
0.0001 



0.0002 
0.0006 
0.0006 
0.0001 
0.0000 
0.0004 
0.0000 
0.0019 
0.0007 



0.0002 \ 
0.0005 
0.0007 
0.0001 
0.0002 
0.0005 
0.0001 
0.0007 
0.0029 j 



The matrix T is a diagonal matrix with diagonal entries given by 

d = (-116.4889, -33.4063, -11.8389, -21.1723, -11.9582,-161.2178, -34.7884, -19.7664, -27.2993) 

and 9 = —31.2689. The eigenvalues of the H> matrix are given by the vector 

e = (-.05025, -.1456, -.0605, -.0424, -.0231, -.0.0053, -.0101, -.0154, -.0131). 

The following Table provides some R and V values that were found as numerical solutions to the equation 
G(R) = a, for selected a's. 



a 


0.05 


0.025 


0.01 


R 


0.9160 


1.0038 


1.1128 


V 


31.6883 


31.7727 


31.8881 



5.2 Example with n > and n + > 

We consider a portfolio that contains call options and put options on equities, so that price of the portfolio 
at time t, is given by: 

5 10 

n(t) = - SiSi{t)] -J2 p At, Sj(t)) (Sj i)s,.(t)] 

i=l j=6 

The prices of each of the options are taken from data from the French CAC 40 Market, and are given in the 
following table: 



Table 2: Data for Ten CAC 40 European Call and Put Options 





Exercize Price 


Interest Rate 


Maturity 


Underlying Price 


Call-BNPPARIBAS 


30.00 


0.05 


3 months 


39.75 


Call-BOUYGUES 


19.00 


0.05 


3 months 


27.30 


Call-CAP GEMINI 


20.00 


0.05 


3 months 


24.00 


Call-CREDIT AGRICOLE 


10.50 


0.05 


3 months 


14.80 


Call-DEXIA 


9.00 


0.05 


3 months 


9.38 


Call-LOREALL 


40.00 


0.05 


3 months 


62.90 


Put-SOCIETEGENERALE 


50 


0.05 


3 months 


64.00 


Put-TFl 


18.00 


0.05 


3 months 


22.02 


Put-THOMSON 


9.00 


0.05 


3 months 


17.13 


Put- VIVENDI 


9.00 


0.05 


3 months 


17.00 
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In this example, using the three month historical data for the ten CAC 40 equities, with the exponential 
moving weighted average (EMWA) and A = 0.94, we obtained the following £. 



S = 



/ 


0.0016 


-0.0001 


0.0012 


0.0005 





0008 


0.0007 


-0.0000 


0.0008 


0.0002 





0002 \ 




-0.0001 


0.0008 


0.0004 


-0.0001 





0000 


0.0000 


0.0004 


-0.0000 


0.0005 





0005 




0.0011 


0.0004 


0.0035 


0.0006 





0010 


0.0007 


0.0004 


0.0013 


0.0006 





0006 




0.0005 


-0.0001 


0.0006 


0.0005 





0002 


0.0004 


-0.0001 


0.0003 


0.0001 





0001 




0.0008 


0.0000 


0.0010 


0.0002 





0015 


0.0007 


-0.0000 


0.0008 


0.0000 





0002 




0.0007 


0.0000 


0.0007 


0.0004 





0007 


0.0010 


-0.0002 


0.0005 


0.0003 





0004 




-0.0000 


0.0004 


0.0004 


-0.0001 


-0 


0000 


-0.0002 


0.0015 


-0.0002 


0.0008 


-0 


0001 




0.0008 


-0.0000 


0.0013 


0.0003 





0008 


0.0005 


-0.0002 


0.0012 


0.0000 





0001 




0.0002 


0.0005 


0.0006 


0.0001 





0000 


0.0003 


0.0008 


0.0000 


0.0018 





0007 


\ 


0.0002 


0.0005 


0.0006 


0.0001 





0002 


0.0004 


-0.0001 


0.0001 


0.0007 





0026 / 



The matrix T is a diagonal matrix with diagonal 

d = (24.186, 8.6269, 21.7320, 4.3111, 15.4949, -4.5815, -82.2915, -22.2079, -1.2957, -1.2822), 

and 9 = —3.8596. The eigenvalues of the D matrix are given by the vector 

e = (-0.1251, -0.0115, -0.0030, -0.0014, -0.0006, 0.0014, 0.0092, 0.0124, 0.0290, 0.1271). 

The following Table provides some R and V values that were found as numerical solutions to the equation 
G(R) — a, for selected a's. 



a 


0.05 


0.025 


0.01 


R 


0.6069 


0.7176 


0.8455 


V 


4.0438 


4.1171 


4.2166 



6 Conclusions 

We have considered an application in the domain of Risk Management to estimate the Value-at-Risk of 
the portfolio of an option. We have reduced the problem to a multiple integral over a hyperboloid. This 
type of integral can be approximated using techniques described by Genz and Monahan [7] and Sheil and 
O'Muircheartaigh |11| . 
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